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Abstract 

An outstanding challenge for models of non-perturbative quantum gravity is the consistent 
formulation and quantitative evaluation of physical phenomena in a regime where geome- 
try and matter are strongly coupled. After developing appropriate technical tools, one is 
interested in measuring and classifying how the quantum fluctuations of geometry alter the 
behaviour of matter, compared with that on a fixed background geometry. 

In the simplified context of two dimensions, we show how a method invented to analyze 
the critical behaviour of spin systems on flat lattices can be adapted to the fluctuating ensem- 
ble of curved spacetimes underlying the Causal Dynamical Triangulations (CDT) approach 
to quantum gravity. We develop a systematic counting of embedded graphs to evaluate the 
thermodynamic functions of the gravity-matter models in a high- and low-temperature ex- 
pansion. For the case of the Ising model, we compute the series expansions for the magnetic 
susceptibility on CDT lattices and their duals up to orders 6 and 12, and analyze them by 
ratio method, Dlog Pade and differential approximants. Apart from providing evidence for 
a simplification of the model's analytic structure due to the dynamical nature of the geome- 
try, the technique introduced can shed further light on criteria a la Harris and Luck for the 
influence of random geometry on the critical properties of matter systems. 



^email: d.benedetti@phys.uu.nl, r.loll@phys.uu.nl 
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1 Coupling quantum gravity to matter 



If a quantum theory of gravity is to describe properties of the real world, it must tell us if and 
how matter and spacetime interact at extremely high energies. Current, incomplete models 
of quantum gravity usually approach the problem by first trying to construct a quantum 
theory of spacetime's geometric degrees of freedom alone, and then adding matter degrees 
of freedom in some way. One is then interested in whether and how the "pure" gravity 
theory gets modified and how the dynamic aspects of quantum geometry may influence the 
matter behaviour. A priori, different scenarios at very short distance scales are thinkable: 
the behaviour of quantum geometry may be completely dominant, geometric and matter 
degrees of freedom may become indistinguishable, or the matter-coupled theory may have 
no resemblance with the pure theory at all. 

In practice, our knowledge of non-perturbative models of quantum gravity, with or with- 
out matter, and whether they possess physically interesting properties, is at this stage lim- 
ited. However, the development of some of the purely geometric models has advanced to 
a point where one can reasonably consider their more complex matter-coupled versions, in 
order to obtain a first quantitative idea of their behaviour. Such attempts have to over- 
come a number of serious difficulties. First, there is the problem of what to calculate in a 
background-independent formulation, that is, to come up with geometric invariants which 
relate to physical observables. Second, there is the problem of how to calculate, once such 
quantities have been identified, that is, the need to establish a viable calculational scheme. 
In case one uses a numerical approximation, say, to evaluate a non-perturbative path in- 
tegral over geometries by Monte Carlo methods, one has to deal with the computational 
limitations of the hardware. In case one uses some expansion scheme, it must have non- 
standard features in order to be inequivalent to the usual (non-renormalizable) gravitational 
perturbation series. 

Many difficulties stem from the fact that in the case of quantum gravity, conventional 
calculational methods have to be adapted to a situation where there is no fixed background 
spacetime, and instead "geometry" is among the dynamical degrees of freedom. In this 
paper we will look in detail at an example of such a generalized method for a particular class 
of quantum systems of gravity and matter. We will work in the simplified context of two 
spacetime dimensions, where it is fairly clear what one would like to calculate, and where 
the computational limitations are minimal compared with the full, four-dimensional theory. 
More specifically, we will look at spin systems coupled to dynamical Lorentzian geometries, 
given in the form of an ensemble of triangulated spacetimes in the framework of Causal 
Dynamical Triangulations (CDT) (see [1] for recent reviews). 

We will show how a time-honoured method of estimating the critical behaviour of a 
lattice spin system, namely, by considering finitely many terms of a weak-coupling (i.e. a 
high-temperature) expansion of its thermodynamic functions can be adapted to the case of 
"quantum-gravitating" lattices. The method consists in classifying local spin configurations 
(which can be represented by diagrams consisting of edges or dual edges of a triangulated 
spacetime), and determining their weight, i.e. the probability to find such a configuration in 
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the ensemble of dynamical triangulations. Diagrammatic techniques have also been employed 
recently in another coupled model of geometry and matter in three dimensions, in the context 
of a so-called group field theory, an attempt to generalize matrix model methods to describe 
ensembles of geometries in dimension larger than 2 (see [2] for a recent review). 

The model we will be studying in detail is that of an Ising model whose spins live on 
either the vertices or triangles of a two-dimensional simplicial lattice and interact with their 
nearest neighbours. Summing over all triangulated lattices in the usual CDT ensemble, each 
with an Ising model on it, gives rise to the matter-coupled system whose properties we are 
trying to explore. An exact solution to this Lorentzian model has not yet been found. This 
has to do with the fact that in terms of the randomness of the underlying geometry the model 
lies in between that on a fixed regular lattice and that on purely Euclidean triangulations, 
but is sufficiently different from either to make the exact solution methods known for these 
cases inapplicable. 

What underlies the successful application of powerful matrix model methods to solve 
the model of Euclidean dynamical triangulations with Ising spins [3j is the fact that on 
Euclidean geometries, no directions are distinguished locally, which matches with the fact 
that the matrix model generates configurations (in the form of triangulated surfaces or their 
dual graphs) which are unconstrained gluings of its elementary building blocks. However, 
it turns out that for the purposes of quantum gravity this property of local "isotropy" is 
simply not good enough. Ensembles of unconstrained gluings of triangular building blocks 
(so-called four-simplices) in four dimensions, even if they are required to be manifolds of a 
fixed topology, have been found to be dominated by highly degenerate polymeric or tightly 
clustered geometries, neither of them making good candidates for a theory of quantum 
gravity (see [1] for a review). This was precisely the reason for introducing the Lorentzian 
CDT version of the original (Euclidean) dynamical triangulations, which does keep track of 
the non-isotropic light-cone structure and fiow of time on its spacetimes. Remarkably, in this 
model, evidence for the desired four-dimensionality of its geometric ground state on large 
scales has recently been found [3, El E] . In more recent times, the introduction of somewhat 
different "causality" constraints has been advocated in a Euclidean spin-foam approach to 
quantum gravity (see [8] and references therein). 

Returning to our previous argument, it is notoriously difficult to introduce additional 
restrictions concerning the local nature of geometry into a matrix model without destroying 
its simplicity and, more importantly, its solubility. This also seems to apply to the two- 
dimensional CDT model, whose pure-gravity version was first solved exactly in M by transfer- 
matrix method, with the Ising-coupled model studied subsequently in [101 E] O In the latter 
case there is strong evidence from numerical simulations that the spin model behaves in the 
presence of gravity (at least as far as its critical exponents are concerned) just like it does on a 
fixed, regular lattice. This is supported by a "semi-analytical" high-temperature expansion, 
whose initial results were merely quoted in the original work of [TU]. In the present paper, we 

^We would expect a similar issue to appear in the group field theory models mentioned earlier, which 
beyond triangulated manifolds generate a vast number of geometric configurations, many of which one 
presumably would like to get rid of by imposing suitable constraints (see also the remarks in |2]). 
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give a detailed technical exposition of how to perform the relevant graph counting on causal 
dynamical triangulations, both for the standard Ising model and its dual. The algorithms 
found have allowed us to compute the diagrammatic contributions to the Ising susceptibility 
to order 6 and order 12 for CDT lattices and their duals, respectively. 

The pure counting results were announced in a recent letter p2], where we put forward 
the conjecture that the best way to extract information about the universal properties of a 
matter system on flat space is by coupling it to an ensemble of quantum-fluctuating causal 
triangulations! The strongest evidence for this conjecture so far comes from the fact that 
an evaluation by straightforward ratio method of a rather small number of susceptibility 
coefficients already leads to results in surprisingly close agreement with the exactly known 
value for the susceptibility. As we will also reiterate in the present paper, the evidence for 
this conjecture is still too limited to reach any definite conclusion. It is an issue that is being 
studied further. 

Generally speaking, one can look at our method and results also from a purely condensed 
matter point of view, as an example of how the universal properties of spin and matter 
systems are affected by introducing random elements into its definition. For example, the 
inclusion of random impurities in a lattice model leads to Fisher renormalization [13] of 
the critical exponents in the so-called annealed case, where the disorder forms part of the 
dynamics. By contrast, in the quenched case the Harris-Luck criterion is conjectured to hold 
[Ti] . To what extent these ideas extend to the case of random connectivity of the lattice 
- which is the one relevant to non-perturbative quantum gravity models - is still not very 
clear, despite a number of results in this direction, involving spin models with (quenched) 
geometric disorder coming from Euclidean dynamical triangulations and Voronoi-Delaunay 
lattices (see, for example, fTEl \TE\ ). 

One of the reference points for such investigations is that of the two-dimensional Ising 
model on fixed, regular lattices, which in absence of an external magnetic field can be solved 
exactly in a variety of ways (see [13, [H]). In this case, the critical exponents, which charac- 
terize the system's behaviour near its critical temperature, are given by the so-called Onsager 
values a = 0, (3 = 0.125 and 7 = 1.75 for the specific heat, spontaneous magnetization and 
susceptibility. In sharp contrast, in the case of the Ising model on Euclidean dynamical 
triangulations mentioned earlier the disordering effect is so strong that the same critical 
exponents are altered to a = —1, P = 0.5 and 7 = 2 [3J. On the other hand, despite the 
annealed geometric randomness present, the corresponding Ising model coupled to causal 
dynamical triangulations seems to share the flat-space exponents, indicating a perhaps sur- 
prising robustness of the Onsager universality classH 

In the present piece of work, we will concentrate for simplicity and deflniteness on the 
CDT-Ising system, although our method should also apply straightforwardly to other spin 
models and may inspire similar diagrammatic expansions in other, discrete models of quan- 
tum gravity. In the next section we will review the main features of causal dynamical 
triangulations in two dimensions. Sec. 3 introduces the high-temperature expansion on reg- 

■^The study of fermions coupled to two-dimensional causal dynamical triangulations was initiated in [19j . 
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ular and dynamical lattices. In Sec. 4, after recalling some general definitions and results of 
graph theory, we present a set of rules which in principle allow us to compute the weight of 
any susceptibility graph on the fluctuating CDT lattice and its dual, and explain the method 
with the help of illustrative examples. The finite number of terms of the series expansions 
we obtain in this way are analyzed in Sec. 5 with the help of the ratio, Dlog Pade, and 
differential approximants methods. Sec. 6 contains some remarks on the extension of our 
method to low-temperature expansions, before we summarize our findings in Sec. 7. 

2 Causal Dynamical Triangulations revisited 

Despite the apparent failure of describing gravity as a renormalizable quantum field theory, 
the hope that it may nevertheless admit a well-defined non-perturbative formulation has been 
both inspiration and justification for several approaches to quantum gravity, like Dynamical 
Triangulations [20], Loop Quantum Gravity [211 122], Exact Renormalization Group [23], 
Spin Foams ^24] and Causal Sets |25j . 

The way to a non-perturbative quantization pursued by Dynamical Triangulation is that 
of giving meaning to the formal (EuchdeaiJJ path integral 



for the gravitational action, consisting of the usual Einstein-Hilbert term, a cosmological 
term, and possibly others. To accomplish this task one starts from an initially discrete 
formulation where the sum over geometries, i.e. the integral over metrics modulo diffeomor- 
phisms, is defined as the continuum limit of a sum over simplicial manifolds T. Roughly 
speaking, a simplicial manifold is a collection of d-dimensional triangular building blocks 
called "simplices" (generalizing triangles in d = 2 and tetrahedra in d = 3), glued together 
along their {d — l)-dimensional faces (or edges in ci = 2) in such a way that local neighbour- 
hoods look (i-dimensional. To each simplicial manifold T one can associate a Boltzmann 
weight e""^-^, where the Regge action Sii{T) represents a discrete version of the Euclidean 
gravity action. The simplicial building blocks are taken to be all identical, with the length 
of each one- dimensional edge fixed to a cut-off value a. The thus regularized path integral 
takes the form 



where C{T) is a symmetry factor. Note that the diffeomorphism invariance of the contin- 
uum formulation has been taken care of by adopting an explicitly coordinate-independent 
formulation. (We refer to [2S] for details of motivation and construction.) 

"^An important feature of the causal version we will be using is that one starts with Lorentzian signature 
and a complex path integral Z — / 2?[.gpi/]e*'^^"[^f"'l, and then performs a Wick rotation to the Euclidean 
version ([1]), which still retains a memory of the causal light-cone structure of the original spacetimes. For 
the purposes of the current presentation we will work with the Euclideanized version of the Lorentzian CDT 
path integral throughout. 






(2) 
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This approach requires two restrictions on the ensemble of simphcial manifolds to make 
it a viable candidate for a four-dimensional theory of quantum gravity. First, the spacetime 
topology must be fixed, since otherwise the number of inequivalent geometries grows facto- 
rially with the volume and the Boltzmann weight would have no chance of controlling the 
resulting divergence. The second restriction has emerged only after a detailed study of the 
dynamics of the dynamically triangulated model. It turned out that the original, purely Eu- 
clidean model behaved too wildly to admit a sensible classical limit. In particular, no region 
in the phase space of the underlying statistical model of random geometry was found where 
the large-scale geometry was extended and four-dimensional (as, for example, measured by 
its large-scale Hausdorff dimension dn)- 

The restriction proposed in [9] to resolve this impasse was to impose a causal structure on 
the triangulated spacetimes in such a way that the simplices interpolate between successive 
spatial sliced of constant integer proper time. The ensemble of geometries with these two 
restrictions implemented underlies the approach of Causal Dynamical Triangulations. So 
far, CDT have lived up to their original expectation: to the best of current knowledge, 
the quantum geometries generated dynamically from simplicial building blocks of dimension 
d = 2,3 and 4 have an effective large-scale dimension equal to (i [HI [23, El EB]- 

In this paper we shall restrict ourselves to the case d = 2, which brings about several 
simplifications for the CDT model, including a graphical one, in the sense that triangulated 
spacetimes from the causal ensemble just introduced can be represented by planar drawings, 
like the one in Fig. [H In the same figure we also show the corresponding dual graph, where 
triangles have been substituted by trivalent vertices and the gluing of neighbouring triangles 
by links connecting the corresponding vertices. The representation by dual graphs has been 
used in some analytical investigations [25], and will be used in the present article when we 
discuss the "dual" Ising model with the spins placed on triangles instead of vertices. 

A second and more important simplification in the two-dimensional case is that the 
Einstein-Hilbert term (the integrated scalar curvature) is purely topological and can therefore 
be dropped from the action, leaving only the cosmological term proportional to the volume 
of the manifold. This means that in our approach in d = 2 the path integral ([T]) is replaced 
by the partition function 

N Tn 

where labels distinct (causal) triangulations of fixed volume A^, the number of triangle^, 
and A is the bare cosmological constant. Asymptotically can be shown to behave like 

Z = ^e(^-^)^+°(^), (4) 

N 



slice is defined as a (d — l)-dimcnsional spatial submanifold of constant topology S of the simplicial 
spacetime manifold. This fixes the spacetime topology to be of type S x M. For definiteness, we will use 
S = 5^ in what follows. 

^With all triangles equilateral of side a, the volume of each triangle is -^a^, and hence the total volume 
proportional to N. 
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Figure 1: Example of a triangulated piece of spaeetime in the CDT model (left) and the corresponding dual 
graph (right). The picture captures the way triangles are glued together, but does not represent faithfully 
intrinsic distances, and therefore the curvature properties of the simplicial geometry. If we were to respect 
this property we would not be able to draw the picture in a plane, because the model's triangles are all equal 
and equilateral. 

with Ac = log 2. It is straightforward to see that the model is well defined for A > Ac and 
that in the limit of A ^ Ac the average discrete volume (N) goes to infinity. In this limit 
one can send the cut-off a to zero while keeping the physical volume V — ^a^N finite, thus 
obtaining the continuum limit of the model. 

3 The Ising model and its high-T expansion 

Given any lattice G of volume N with v vertices and / links (edges), the Ising model on G 
in the presence of an external magnetic field h is defined by the partition function 



where is the Hamiltonian of the spin system, /3 l/{kBT) the inverse temperature, 

(Tj the spin variable at vertex taking values 0"^ = ±1, and {ij) denotes nearest neighbours. 
In standard notation, we will use K = j3J, where J > is the ferromagnetic spin coupling, 
and H = jSmh, where ni is the spin magnetic moment. The fact that a^aj = ±1 allows us 
to write 





(5) 



{<Tf=±l}igG 



{a-i=±l}fgG 




{1 + u cTiaj) cosh(i^'), 

(l+TCTj)cOsh(ii"), 



(6) 
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in terms of newly defined variables u = tanh(i^) and r = tanh(if), and rewrite the partition 
function as 

= cosli'(if ) cosli''(/7) J2 Ili^ + « ^i^j) Ili^ + ^ 

= cosh^(i^') cosh''(if) ^[{)- + a^aj + ^ o-iCTjcrfca/ + ...) x (7) 

Wi} (ij) {i3),{kl)d^l 

X {l + T^ai + T^^aia^ + ...)] 

= cosh'(K) cosh''(i/) r). 

Because of the sum over spin values ±1 it is clear that every term containing a spin at a 
given vertex to some odd power will give a vanishing contribution. Thus we need to keep 
only terms with even powers of cTj's, which are terms where each ai belongs either to an even 
number of nearest neighbour couples (TiCTj (from the u-terms) or to an odd number of them 
and to one of the spins coming from the r-terms. It is not difficult to convince oneself that 
each such term is in one-to-one correspondence with a graph drawn on the lattice whose 
length is given by the power of u and whose number of vertices with odd valence is given by 
the power of r. The function F^^\u, r) introduced in eq. ([7]) is then the generating function 
for the number of such graphs that can be drawn on the lattice, and the factor 2^ comes from 
the sum over the spin configurations. The representation we have obtained in this manner 
is a high-temperature expansion, since for infinite temperature T we have u = r = 0. We 
will return to this graphic interpretation in the next section, after having recalled some basic 
definitions from graph theory. In what follows, we will only be interested in the case of 
vanishing external magnetic field. 

All we have said up to now does not require any specific properties of the lattice G, but 
works for any lattice, regular or not. This makes it straightforward in the framework of CDT 
to couple the Ising model to gravity. We simply associate a complete Ising model with each 
triangulation of the ensemble, viewed as a lattice, by putting the spins at the vertices of the 
triangles, and then perform the sum over such triangulations, leading to 

Z{K,H,N) = Y^ Yl e-^'^^'^^l. (8) 

The situation where the volume is allowed to fiuctuate is also of interest in a quantum gravity 
context and described by the grand canonical partition function 

ZiK,H) = J2e-'''ZiK,H,N), (9) 

N 

where the role of chemical potential is played by the cosmological constant A. As usual in 
dynamically triangulated models, the asymptotic behaviour as function of of the canonical 
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partition function in the infinite-volume limit N ^ oo will determine the critical value of A 
at which the continuum limit can be performed. 

In trying to understand the behaviour of the Ising model coupled to quantum gravity in 
the form of causal dynamical triangulation, we will make crucial use of the known probability 
distribution P{q) of the coordination number q of vertices (the number q of links meeting 
at a vertex) for the pure gravity case. This is relevant because we will be expanding about 
the point = 0, at which the spin partition function ([7]) reduces to a trivial term 2^ and 
the geometry is therefore that of the pure 2d Lorentzian gravity model. The probability 
distribution in the thermodynamic limit was derived in [TO], resulting in 

^(9) = 1^- (10) 

The main idea is to use this distribution, together with some information about correlations 




Figure 2: Comparing the probability distributions P{q) of the vertex coordination numbers q for three 
different types of random triangulations: Euclidean dynamical triangulations [30^, Poissonian Delaunay 
triangulations fWD and causal dynamical triangulations [TU] , 

between distributions at different vertices, to re-express the sum over triangulations in the 
high-temperature expansion of the partition function, ([H]), by an average over coordination 
numbers. In terms of geometric randomness, the pure CDT model lies in between the 
Voronoi-Delaunay triangulations based on Poissonian random distributions of vertices and 
the planar triangulations underlying the approach of Euclidean dynamical triangulations. A 
comparative plot of the probability distributions for the coordination numbers in the three 
different types of geometry is shown in Fig. [21 



9 



Let us point out a further subtlety which arises from the fact that the underlying lattices 
are not fixed, but fluctuating, and therefore do not have a fixed "shape" . In this case it can 
happen that in the graph counting for the high-temperature expansion of some thermody- 
namic quantity topologically non-trivial graphs must be taken into account. By this we mean 
closed non-contractible graphs which wind around the spacetime. This does not invalidate 
the method in principle, but requires a more detailed knowledge of the global properties of 
the spacetime. In two-dimensional CDT, where the spacetime topology is usually chosen 
to be that of a cylinder (a spatial circle moving in time), this would be closed loop graphs 
which wind one or more times around the spatial S^, and might be as short as a single link. 
Such pinching configurations certainly exist, but their contribution to the graph counting 
has been shown to be irrelevant in CDT, because it is subleading in in the thermodynamic 
limit [lU]. 

4 Graph embeddings 

4.1 Terminology of graph theory and embeddings 

To prepare the ground for our counting prescription, we will first review some of the relevant 
terminology of graph theory and embeddings, along the lines of reference [32] . 

A linear graph is a collection of v vertices and / lines (or edges) connecting pairs of 
vertices. A simple graph is a graph in which two vertices are connected by at most one 
line and in which single lines are not allowed to form closed loops. A graph is said to be 
connected if there is at least one path of lines between any two vertices, and disconnected if 
for some vertex pair there is no such path. 

The cyclomatic number c of a connected graph g is defined as 



and represents the number of independent cycles in the graph. 

The degree (or valence or coordination number) of a vertex is the number of edges incident 
on that vertex. 

Two graphs are said to be isomorphic if they can be put into one-to-one correspondence 
in such a way that their vertices and edges correspond. They are called homeomorphic if 
they are isomorphic after insertion or suppression of any number of vertices of degree 2 (this 
operation being defined on an edge in a trivial way). 

Homeomorphs are thus graphs with the same topology, in particular, with the same 
cyclomatic number. It is then possible to classify graphs in terms of irreducible graphs. 
An irreducible graph is one in which all vertices of degree 2 have been suppressed. (Note 
that in general it will no longer be a simple graph since it will contain single-line loops.) 
Typical irreducible graphs whose homeomorphs we will encounter are tadpoles, dumbbells, 
figure eights and 6 -graphs (Fig. [3]). 

In the high-temperature expansion of an Ising model it is useful to think of the lattice 
with its V vertices and / edges as a graph G. The graphs appearing in the expansion of the 



c{g) = l-v + l, 
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(a) (b) (c) (d) 



Figure 3: Four kinds of homeomorphs: (a) tadpole, (b) dumbbell, (c) figure eight, and (d) 6'-graph. 



thermodynamic function of interest are then graphs embedded in G. A general definition of 
embeddings of graphs is the following. Let G be a general graph. A graph if is a subgraph 
of G when all its vertices and edges are vertices and edges of G. An emheddin^ of a graph 
(yf in G is a subgraph G' of G which is isomorphic to g. 

The lattice constant (g;GY^^ of a graph g on G is defined as the number of different 
embeddings of g in G. In practice it is often useful to work with the lattice constant per site, 
defined as {g;G) = ^{g;G)'^^l 

With these definitions in hand we can now rewrite the function F^^^ (u, r) introduced in 
(CD as 

oo 

F(-)Kr)^l + 5:5:<inV-, (12) 

s=0 l>0 

where D^^^ are the lattice constants of graphs of length / and with 2s odd vertices, for which 
F^^\u,t) serves as a generating function. It follows from the extensive nature of the free 
energy of the system that in the thermodynamic limit 

FW(u, t) e^®("'^) = 1 + NOiu, t) + OiN^), (13) 

where the function Q{u,t) does not depend on N. From the usual combinatorial theory of 
non-embedded graphs it is well known that the logarithm of the generating function is the 
generating function for connected graphs. Because of the embedded nature of the graphs 
in the case at hand, the disconnected graphs will still contribute to log F^^^ (-u, r) with a 

"repulsion" term of sign (^_l)(#con"ecied components)~l 

In our analysis of the Ising model on CDT we will be looking at a quantum-gravitational 
average of the magnetic susceptibility at zero external field. For the usual Ising model the 
susceptibility x per unit volume is given by 



1 dHnZ^ 1 , 2 EiD 



L2 



U 



^^"^ N dm |h=o 2^ Ni + Y.MVu'^ ^^^^ 



^To be precise, this defines a weak embedding. By contrast, a strong embedding of g in G is defined as 
any section graph of G which is isomorphic to g, where a section graph is a subgraph of G consisting of 
a subset A of vertices and all the edges which connect pairs of (nearest-neighbour) vertices of A. Since all 
our calculations will concern weak embeddings and the corresponding weak lattice constants, we will from 
now on drop this specification. 
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where in the second step we have substituted in the high-temperature expansion of the 
partition function. By virtue of (fT^ it is clear that the susceptibihty is independent of 
N in the infinite- volume limit (thus justifying our notation x{i^))y ^ind that moreover the 
denominator in the last term in ([T3j) does not contribute at lowest order in N, which is the 
one relevant to the computation. To calculate the susceptibility, it is therefore sufficient to 
compute the term of order N in Yli ^fl^'^^ is well known. 

The computation in the gravity-coupled case proceeds completely analogously. Perform- 
ing the sum over all triangulations at fixed N , as in eq. ([H]), and then letting N oo, one 
obtains 

2 E,(ET.fl,?(r.))«' 2 E,(A',f)n' 
'"ET„i + E,(ET„<'(r.))^' ^i + E,«V 

for the susceptibility in presence of the gravitational CDT ensemble, where we have dropped 
the irrelevant constant term and introduced an obvious notation for the ensemble average in 
the second step. Whereas it was fairly straightforward to verify the cancellations of terms 
of higher order in N between numerator and denominator in the pure Ising case in eq. ( fT4l) 
by an explicit calculation, the analogous computation does not seem feasible in the gravity- 
coupled case, although we know from the same general arguments that it must be realized 
here too. Again the terms in the denominator of the expressions on the right-hand side 
of (fT5l) do not contribute at lowest order, which - at least for this particular observable - 
implies that its annealed and quenched gravitational averages coincide. 

4.2 Useful results on lattice constants 

We now want to identify which graph topologies appear in the high-temperature expansion 
of the magnetic susceptibility. For convenience we can divide our set of graphs into two 
subsets, graphs with zero cyclomatic number (also called Cayley trees), and graphs with 
c > 0. Since the graphs we encounter in the susceptibility series must have two and only two 
odd vertices, it follows that the first subset will only contain connected non-selfintersecting 
open chains (in other words, self-avoiding walks), while the second will contain everything 
else. The connected graphs in the second subset are graphs with loops and one or no open 
ends, while the disconnected ones take the form of a union of one of the previous two kinds 
with any number of even graphs (graphs in which every vertex has even degree). 

The problem we meet in the high-temperature expansion is the evaluation of the lattice 
constants of such graphs. It turns out that there are two kinds of theorems which can 
be applied profitably in the context of our random triangulated lattices. The first is very 
general and reduces the calculation of the lattice constants for disconnected graphs into that 
for connected graphs. This makes it possible to compute the lattice constants for the second 
subset (c > 0) in a unified way. The second theorem is less general but powerful, since it 
gives us a recurrence relation for the expansion coefficients of the susceptibility. 



Xcdt(^^) 
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The reduction theorem for disconnected graphs states that if gi and gj are two graphs 
gi 7^ gj and G is any graph, then 

{gi U gj; G) = {gf, G){gf G) - ^{gi + gj = gk]{9k] G), (16) 

k 

where gi U gj stands for the disjoint union of the two graphs, {gi + gj = gk} is the number 
of possible choices of embeddings of gi and gj in gk having gk as their sum graph, and the 
summation is over all graphs gk obtainable in this way. If gi = gj we can compute the right- 
hand side of (fTBI) as if the two graphs had different colours and thus obtain 2{gi U gi]G). 
If gi or gj are themselves non-connected graphs we can iterate the theorem. In this way 
the lattice constant of the disconnected union of n connected graphs can be expressed as a 
polynomial of order n in the lattice constants of connected graphs (see Fig. H] for some simple 
examples) . 





Figure 4: Two examples of the reduction theorem. The graph G has been dropped in the notation. 

The counting theorem for the susceptibility series was first given by Sykes [33] in 1961, but 
only later proved by Nagle and Temperley [31] as a special case of a more general theorem 
of graph combinatorics. Closer examination reveals that the proof makes no reference to the 
structure of the lattice, but only to its coordination number g, which will turn out to be 
useful when considering the lattice dual to the triangulation. The theorem states that the 
high-temperature series for the susceptibility can be written as 

X{u) = 1 + (1 - au)-^ [quil - au) - 2(1 - u^)S^ + 8(1 + u)\S2 + S^)] , (17) 
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where g = o" + 1 is the coordination number, and the functions Si are defined as 



5i = ^ IdiU^ 



l>0 



l>0 



'^wdi{w) 



u 



«>0 



rr'mi(r, r' 

,r<r' 



The coefficient di is the sum of lattice constants of even graphs with I fines, wfiile di{v) is 
restricted to even grapfis witfi / fines and vertices of degree 2qi, 2q2, characterized by the 
weight w = Yli^ii^i ~ l)/2- The function mi{r,r') is the sum of lattice constants of graphs 
with exactly two odd vertices of degrees 2r + 1 and 2r' + 1. Substituting the left-hand side 
of formula (fTTl) by the expansion = 1 + ^„>o ^nU"^, one finds a recursive relation for 
tfie coefficients a„, namely. 



ai - 2oai_i + o'^ai_2 = 2{l - 2)di-2 - 2ldi+ 



+ 



+ 



wdi{w) + 2 ^ wdi_i{w) + ^ wdi_2{w) 



u)>0 



+ 



(19) 



rr'mi{r, r') + 2 rr'mi_i{r, r') + rr'm;_2(r, r') 



r<ir' 



As a result, at every new order / in our calculation the only new lattice constants we have to 
evaluate are those for graphs with no vertex of degree 1, reducing the computational effort 
considerably. 

This formula can in principle be applied to any non-regular lattice whose vertices are all 
of the same degree. For a fixed non-regular lattice this observation would usually not be of 
much help, because one would still need to compute at each order the new lattice constants 
depending on the complete, detailed information of the lattice geometry. In our case this 
difficulty is not present, since we are summing over triangulations and can simply substitute 
in ( JT9l) the lattice constants averaged over the triangulations, without the need of keeping 
track of the lattice geometry for individual lattices. 



4.3 Counting graphs on the CDT triangular lattice 

Next, we will analyze the evaluation procedure for lattice constants in the quantum-gravita- 
tional CDT model. Computations here are made difficult by the randomness of the coordina- 
tion number. At the outset it is not even obvious that any of the known recurrence relations 
can be used. Our task is to list all possible diagrams, and for each count the number of 
ways it can be embedded in the lattice. On a regular lattice such an operation is tedious 
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but straightforward: starting from a fixed vertex, we trace out all possible sequences of links 
of a given length, say, which do not self-intersect. The discrete symmetries of the regular 
lattice simplify this task greatly. In the simplest case, the counting has to be done only for a 
single initial vertex, with all other graphs obtained subsequently by translational symmetry. 
The analogous counting on the dynamical CDT lattices is complicated by the fact that the 
lattice neighbourhoods of different vertices will in general look different. We will deal with 
this difficulty by setting up an algorithm to count the embedding of a given diagram on the 
ensemble of CDT lattices, making use of the known probability distribution of the vertex 
coordination numbers. 

Before presenting this algorithm we need some more notation. CDT lattices have two 
kinds of links, time-like and space-lik^. In keeping with the usual representation of two- 
dimensional CDT lattices (c.f. Fig. [H left-hand side), we will draw time- and space-like links 
in embeddings of graphs as (diagonally) upward-pointing and horizontal lines, respectively 
(see Fig. [5]). Furthermore we need to keep track of the relative up-down or right-left orienta- 
tion of consecutive links, to distinguish between graphs like (c) and (d) in Fig. [5] (which give 
different contributions to the lattice constant). In addition, the counting for a graph like 
(b) will be identical to that for its mirror images under left-right and up-down reflections. 
These mirrored graphs will be counted by multiplying the embedding constant of the graph 
by an appropriate symmetry factor. A symmetry factor 1 is assigned to graphs symmetric 




(a) (b) (c) (d) 

Figure 5: The possible typologies of embedding, up to reflection symmetry, of a length-2 open chain. The 
angle which the time-like links form with the vertical has no relevance for the graph counting. 

under both up-down and left-right reflections (for example, graphs (a) and (d)), a factor 2 to 
graphs with only one of the two symmetries (like graph (c)) or symmetric with respect to a 
composition of the two (like graph (c) of Fig. [H]), and a factor 4 to graphs without reflection 
symmetry (like graph (6) of Fig. EJor (a) and (6) of Fig. 

The lattice constant {g; T) of a graph g embedded in a graph corresponding to a trian- 
gulation T is computed as 

{g■,T) = J2■^^i9■,T)^, (20) 



^ Even after performing the Wick rotation to Euclidean signature, these two link types remain combi- 
natorially distinguishable because of the special way in which the original Lorentzian simplicial geometries 
were constructed. 
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(a) 



(b) 



(c) 



Figure 6: Three of the different typologies of possible embeddings of a length-3 open chain. 

where {g; T)i is the number of ways a particular typology of embedding {g)i can be realized 
on T, Si is the relevant symmetry factor, and the sum extends over all possible typologies 
of embeddings of the graph g. In the context of quantum gravity, we are interested in the 
average of (!20l) over all triangulations, that is. 



which will take into account the probability distribution of the vertex coordination. 

As mentioned earlier, all calculations will be performed in the thermodynamic limit where 
the cosmological constant is tuned to its critical value Ac = In 2. In this limit, the probability 
of having k incoming time-like links at a particular lattice vertex is given by 



with an identical probability for having k outgoing time-like links at a vertex. Moreover, 
at one and the same vertex the two probabilities are independent of each other. In fact, 
not just for a single vertex are these two probabilities independent, but the same is true 
for all vertices lying in the same space-like, horizontal slice. By contrast, the outgoing 
probability at a vertex will in general condition the incoming probability at another vertex 
on the subsequent horizontal slice. By construction, the probability of having a space-like 
link to the left and to the right of an given vertex is equal to 1. 

Armed with this information we are now ready to compute any of the {g;G)i in fl20p . 
We do not have a general counting formula in closed form, but will formulate a number of 
rules and tools which will enable us to do the counting recursively. We will start with some 
illustrative examples. Examples 3 and 4 will serve as our elementary building blocks in more 
complicated constructions. 

Example 1. The simplest example is that of a length- 1 graph ci, consisting of a single 
horizontal or vertical link. The number of horizontal embeddings is the number of horizontal 
links on a CDT lattice, which by virtue of the Euler relation is given by half the number 
of triangles, i.e. A^/2. If we give the lattice constant per vertex, as is customary in the 
literature, it is therefore {ci)h = 1- By a similar topological argument we can immediately 
deduce the number 2 for the lattice constant of vertical embeddings, but it is instructive to 



(21) 



p{k) 




(22) 
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compute it from the probability distribution (122|) instead. If at a vertex we have k outgoing 
vertical links, there are precisely k ways to embed the one-link graph such that it emanates 
in upward direction from this vertex. Since every vertical embedding is accounted for in this 
way, by assigning the relevant probability and summing over k we obtain 

(^i)^ = E ^ = 2> (23) 

k=l 

in agreement with the earlier argument. The total contribution to the lattice constant at 
order 1 is therefore {ci)h + {ci)^ = 3. 

Example 2. Next, let C2 be a length-2 open chain, and consider the embedding typology 
C2,b shown in Fig. M^b). The result follows straightforwardly from example 1: there is one 
horizontal link per vertex, and there are in the ensemble on average two ways of attaching 
an outgoing vertical line to its right vertex, yielding {c2)b = 2. For the embedding of Fig. 
ISt^a) there is nothing to compute; we have (02)0 = 1- The embedding of Fig. ^d) is the 
first configuration we encounter which extends over two triangulated strips. Because of the 
independence of the probabilities (l22i) in successive strips we obtain 

00 , 

(^^)^ = (E = 4- (24) 

k=l 

Example 3. Now consider the embedding C2,c depicted in Fig. M^c). If the vertex on top 
has k incoming links, there are exactly k{k — l)/2 ways to realize the embedding, leading to 

fe>= = E^|^ = 2. (25) 

fc=l 

Note that the same "inverted- V" embedding on a flat triangular lattice T^g would give 
(c2; Treg)c = 1 instead. Putting everything together and multiplying by the appropriate 
symmetry factors we get {c2)a + ^{c2)b + 2(c2)c + {c2)d = 17 as a total contribution at order 
2. 

In view of more complicated graphs, an alternative and more convenient way of doing the 
calculation in example 3 is by focussing on the probability distributions of the vertices lying 
on the lower space-like slice of the triangulated strip (see also Fig. [Tl^a)). The probability of 
having n space-like links in between the two vertices at the ends of the inverted V is given 
by the probability that there are n — 1 vertices in between with precisely one outgoing link 
each, which is l/2"~^. Thus we obtain again 

(C2;G), = 2 5^^ = 4. (26) 

n=l 
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m 





n 



n 



(a) 



(b) 



Figure 7: Two configurations contributing to the calculations of examples 3 and 4. The labels m and n 
count the numbers of space-like links as indicated. 

Example 4. Consider an open chain of length greater than 2 embedded in such a way 
that only the first and last edges lie along time-like links, with a horizontal chain of m 
space-like links in between. Both time-like links are supposed to lie in the same strip of the 
triangulation, as illustrated in Fig. ^b) . We are going to prove that the lattice constant of 
such a configuration is given by 



The result 4 — 2^""* could be obtained easily by (i) considering the number of ways in which 
each of the two vertices at the top left and right can have an incoming link (giving a total 
of 2 X 2 = 4 possibilities) and (ii) subtracting from that the probability for the two time-like 
links touching each other in their lower vertex, which is 2^~"* (the probability that the m — 1 
intervening vertices have only one incoming link each). However, we rather want to do the 
counting in a way that keeps track of the number of links separating the two vertices at the 
bottom, for reasons that will become clear soon. In other words, we want to determine the 
lattice constant for a closed polygon consisting of two vertical edges, m edges along the top 
and n along the bottom. From a combinatorial point of view this can be rephrased as the 
problem of counting the number of different ways in which n triangles and m upside-down 
triangles can be arranged to form a strip, with the well-known binomial result (™^")- To 
get the lattice constant for our random lattice we still have to include a probabilistic factor 
I for each vertical link in the polygon interior, leading to the factors 2^"™""- in formula ( 1271) . 

The strips decomposition. The embedding typology {g)i of any connected graph g ex- 
tends over a well-defined number of strips in the CDT triangulation, and each part of {g)i 
belongs to a definite strip or horizontal slice. Our first step in the counting procedure is 
to decompose the embedded graphs into strip contributions, thus breaking the graph into 
pieces. The example of Fig. 18] will help us explain the procedure. It shows a particular 
embedding of C4, how it is decomposed into two pieces, and how its lattice constant can be 
computed accordingly. All we have to do is multiply the probability of C2,c (example 2) with 




(27) 
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m links at the bottom with the probabihty of Cm+2,b (example 3), and then sum over m, 
resulting in 

°° 1 90 
W. = E^(4-2^"'") = y (28) 

m=l 

We can think of this operation as a product of an (infinite-dimensional) vector and matrix, 





Figure 8: Example of a strip decomposition. The embedding on the left is split into two strips, after which 
the strip contributions for given m are multiplied pairwise and summed over. 



K2{m,n) 



1 



2"!— 1 ' 

1 ( m + n 



2m+n-l \ jji 



(29) 



followed by taking the trace in order to impose the open boundary condition, that is, 

(C4). = Tr(iri-/^2). 



(30) 



Fig. [H] is a simple case with only two vertical lines in each strip. In general there will be 
more lines, like in Fig. However, there are only two possibilities for pairs of neighbouring 
vertical links: they can be either disjoint or have exactly one vertex in common, so that the 
probability associated with them in a particular embedding in a particular triangulation will 
be given by either K2{m, n) or Ki{n). We can then associate a probability to a given pattern 
of vertical lines in a strip with fixed distance between them (in terms of the horizontal links 
in the in- and out-slice), which is the product of the corresponding Ki- or if2-probabilities, 
times a factor | for each shared vertical line between two neighbouring i^j-patterns in the 
same strip (see Fig. [9l). In the end we have to glue the strips back together again, with 
conditions such as to get the desired graph, and sum over the allowed values of horizontal 
lengths. As an illustrative example, we obtain for Fig. the multiple sum 

oo oo oo m,2 — l oo -. 

EE E E E[^2(/i+Ji,mi)-ir2(/2,m2)][Ki(mi)-K2(m2-j2,r2)-iri(j2)], (31) 

jl=l mi=l m2=2 j2=l n=l 
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where we have grouped the factors corresponding to the same strip in square brackets. As 
should be clear from this example, the gluing and the range of the summations for the various 
graph embeddings will have to be taken care of on a case-by-case basis. 




n 



Figure 9: Example of a more complicated strip decomposition where each strip consists of a sequence of 
elements of the kind illustrated in Fig. [71 The sum refers to the total sum in the expression ([3T|) which has 
to be performed after the proper weights have been assigned to each element of the graph. 

Using these techniques we have repeated and independently confirmed the calculation of 
the first five orders given in [TDj and then extended it to order 6. Table [T] gives the lattice 
constants (times 2) per number of vertices (which we will call the susceptibility coefficients) 
of the open chains, the other graphs, and of their grand total. The effort involved in this 
extension is by no means minor. There are 387 order-6 open chains that have to be counted 
individually, with the identifications introduced in this section, and for many of them the 
computation requires a careful use of the strip decomposition and a subsequent evaluation 
of the sums involved. 



n 


open 


closed + disconnected 


total (a„) 


1 


6 





6 


2 


34 





34 


3 


174 





174 


4 


859| 


-12 


8471 


5 


4152| 


-1621 


3990 


6 


19800^1 


-1416i 


18384f 



Table 1: Susceptibility coefficients for the CDT lattice up to order 6. The terminology 'open', 'closed' and 
'disconnected' refers to the topology of the graphs. 
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4.4 Counting graphs on the dual CDT lattice 



On the dual lattice (an example of which is depicted in Fig. [H right-hand side) we can apply 
formula ( JT9l) . which for g = 3 reduces to 



Because of the low coordination number the only even graphs that can appear are connected 
or disconnected polygons, and the only contributions to m/ are from ^-graphs and dumbbells 
with or without disconnected polygonal components. For the disconnected graphs we can use 
the reduction theorem. It should be noted that in the overlap decomposition new topologies 
of graphs appear, for example, graphs with more than two odd vertices, still rendering the 
higher-order computations non-trivial. 

It is not difficult to derive the probability distribution relevant for the dual lattice from 
the one of the original lattice. The probability of having two horizontal links at a vertex, 
one to the right and one to the left, is again 1. The probability for an incoming or outgo- 
ing vertical link at a vertex of a given chain of horizontal links is given by the probability 
of having a triangle or an upside-down triangle in the corresponding strip of the original 
triangulation. Since the two cases are mutually exclusive, the in- and out-probabilities are 
not independent and each takes the value 1/2, which may also be thought of as the weight 
associated with a triangle in the original lattice. With these rules the lattice constant of 
an open chain of length 1 is 3/2: 1 for the horizontal link (associated to its left vertex, 
say, in order to avoid overcounting), plus 1/2 for the vertical link. Alternatively, this easy 
example can be computed by use of the Euler relation. For higher-order graphs it is often 
convenient to go back to the original lattice and count - with the appropriate weights - the 
configurations which can be associated to the dual graph under consideration. 

Example 1. Let us consider again a closed polygon with only two vertical links, with m 
links on the top and n on the bottom, but this time on the dual CDT lattice (see Fig. [TOj) . 
A careful evaluation of the probabilities involved yields 



for its lattice constant, where k is the number of vertical links of the lattice which lie inside 
the polygon, and (■)* denotes the average lattice constant per vertex of the dual graph. The 
binomials in flH^ count the number of possible ways to arrange the k internal links, and the 
remaining probability - expressed as a power of 1/2 - turns out not to depend on k. 

As already mentioned, the in-and out-probabilities at a vertex of the dual lattice are 
not independent, which means that we cannot apply the strip decomposition of the original 
lattice directly, but instead have to proceed more carefully. 



ai - iaai-i + 4a/_2 = 2(/ - 2)di^2 - 2ldi+ 

+ 8 [m;(l, 1) + 2m,_i(l, 1) + m;_2(l, 1)] . 



(32) 




k=0 



(33) 
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Figure 10: An embedding of the polygon considered in example 1. 

Example 2. As an example of this, consider a polygon embedding that extends over two 
strips (Fig.dl]), a type of configuration that appears from order 8 onward. Its lattice constant 
is given by 



mi-l / ^\ m2-l 



fci=0 ^ ^ fc2=0 

mm(fci ,711+712—4) 

i=0 



1 



k 



X 



ni+n2- 2\ 1 

i I 2™-i+™'2+".i+".2+A;i-i+A;2+2 ■ 



(34) 



In eq. ( l34l) . mi and m2 count the links on the top and bottom horizontal lines, and ni and 
n2 count the two sets of contiguous links on the central horizontal line. The numbers of 
internal vertical links in the upper and lower strip are denoted by ki and k2, while i counts 
how many links out of the ki end on one of the two intermediate horizontal lines. The logic 
behind the various counting factors appearing under the sums is very similar to that of the 
previous example. 

We have found the method illustrated by examples 1 and 2 the most compact one to deal 
with the counting of graph embeddings on the dual CDT lattice. It has enabled us to push 
our "counting by hand" to order 12, the results of which are given in Table O Because of 
the fixed, low coordination number of the dual vertices, there are far fewer graphs at a given 
order than there are on the original CDT lattice. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


an 


3 


6 


12 


23 


42+f 


78+i 


142+1 


258 


461+i| 


820+1 


1446+i 


2532+i| 



Table 2: Susceptibility coefficients for the dual CDT lattice up to order 12. 
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Figure 11: An embedding of the polygon considered in example 2. 

5 Series analysis 

5.1 Review of different methods 

Having computed the two high-temperature series to some order, we will now turn to an- 
alyzing them, and try to extract information on the critical properties of the underlying 
gravity-matter systems. The standard analytic methods are well-known from the case of 
regular lattices and reviewed in [3^. We will recall them briefly here in order to be self- 
contained. 

The first and most straightforward tool is that of the ratio method, which works as follows. 
Assuming a simple behaviour of the susceptibility of the form 

x{u)r.A{u)(l--Y\B{u) (35) 

near the critical point Uc, with analytic functions A and B, its series expansion xl"^) = 
„>o (^nu"' should yield (to first order in 1/ n) 

.„.^ = lfl + l^). (36) 



One can then generate sequences of estimates of the critical parameters from sequences of 
point pairs {r„,r„_i} [35], namely, 

n{2 - n)rn + (n - l)V„,_i 

In = 7 TT , (37) 

nrn - [n- l)r„_i 
1 

Uc,n = 7T . (38) 

nr„ - (n - l)r„_i 

In general, these will converge very slowly. In addition, if the series has an oscillatory 
behaviour the method can yield alternating over- and underestimates. In this case a fit of a 
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whole sequence {rmin, ■■■,'>^max}, computed with the help of (!36l) . may be better suited than 
the sequence of estimates to find the straight line masked by the oscillations (r^m is chosen 
properly to exclude large deviations from (IHUI) at small n, and r^ax is the ratio computed 

for n = Umax)- 

Oscillations and other irregularities in the expansion are caused by a more complicated 
behaviour of the thermodynamic function, for example, the presence of other singularities in 
the complex plane close to the physical singularity, which in unfortunate cases may even lie 
closer to the origin than the singularity of interest. In case the behaviour near the physical 
singularity is like 

x{u)r~.A{u)(l--y\ (39) 

with A{u) a meromorphic function with singularities close to Uc, a method that should give 
better results than the ratio method is that of the so-called Pade approximants. The method 
consists in the approximation of a function, known through its series expansion to order jV, 
by a ratio of two polynomials of order C and M.^ subject to the condition C + M. < Af, 

with -Fv(x) = ^^o'^fc^^- common usage, the notation [i3/A^] indicates the order of the 
polynomials used. 

In our specific case one takes as the function F{x) to be approximated the derivative of 
the logarithm of xi^), so that Uc can be recognized as a pole and 7 as the associated residue 
in 

■f logxiu) = ^— (1 + 0{u, - u)) . (41) 
au Uc — u 

This is also referred to as the Dlog Pade method. Usually one only looks at the tridiagonal 
band [(A/" — l)/Af], [Af/Af] and [{J\f + 1)/A/] because of the invariance of the diagonal Pade 
approximants under Euler transformations. 

The Pade approximants work well only when there is no additive term B{u). In presence 
of such a term there is a generalization - known as differential approximants - to account 
for functional behaviours of the form 

n 

x[u)^\[A,(l--y + B{u), (42) 

1=1 

with B{u) and Ai, y4„ analytical functions and Mi, m„ a set of singular points. If we are 
considering the logarithmic derivative of a function f{x) we can rewrite fHUj) as 

Pc{x)f{x) - QM{x)f{x) = 0(x^+^+i). (43) 

The idea of the differential approximants method is to generalize this equation according to 

i?^,(x)F;^(x) + QmA^)Fm{x) - Pc{x)Fj^{x) = S^{x) + 0(x'=+^+-^^+-^^+i), (44) 
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where Rm2{x) and Sic{x) are two other polynomials of order and /C, and where one can 
substitute the differential operator hy D = x-^, forcing the point at the origin to be a regular 
singular point. In the following we will only consider the special case R_M2 ~ ^1 §iving rise 
to the so-called inhomogeneous Ist-order differential approximant, denoted by [K./C] A4i]. 
With this method, the exponent 7 can be evaluated as 

7 = (45) 

where simple roolo of Qmi{x), which gives an estimate of the critical point. 



5.2 The CDT lattice series 

The high-temperature series expansion for the Ising susceptibility on the CDT model is given 
by 

1 16 

Xcbt{u) = 1 + 6u + Mv? + 174m^ + (847 + -)m^ + 3990m^ + (18384 + —)u^ + 0[u^). (46) 

Applying formulas flH7|) and flHHj) we get the sequence of estimates for 7 and reported in 
Table El The results of linearly fitting the whole sequence of ratios {r2, r^ax} instead 



n 


Uc 


7 


3 


0.2489 


1.819 


4 


0.2424 


1.721 


5 


0.2459 


1.791 


6 


0.2438 


1.740 



Table 3: Two-point linear extrapolations of the critical point and the critical exponent of the series 
obtained using ([37|l and p8|l . 



are listed in Table HI and a plot of the ratios {r2, rg} versus 1/n together with their fit 
is shown in Fig. [T21 Despite the shortness of the sequence, the linear fit is of very good 
quality. It is worth noting that if we remove the point r2, which of course is subject to the 
largest deviations from a pure 1/n-behaviour due to higher-order terms, and fit the truncated 
sequence {rs, ...,r6}, we obtain 7 = 1.749 and Uc = 0.2441 (see Table [5]). 

It is surprising how fast the ratio sequence seems to converge to the result 7 = 1.75 of 
the exact solution for the flat regular compared to the analogous one obtained from 

the series expansion on the plane triangular lattice (see [SB] for the expansion coefficients to 
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^max 


Uc 


7 


3 


0.2488 


1.820 


4 


0.2462 


1.789 


5 


0.2458 


1.783 


6 


0.2454 


1.779 



Table 4: Linear fits of the sequences {r2, ■••,?maa;} for tlie series witli assumed functional form ([SH]) . 




0,15 0,20 0,25 0,30 0,35 0,40 0,45 0,50 0,55 



1/n 



Fi gure 12: Plot of the ratios (j36p for the case of the CDT lattice (data from Table H)) . 



^max 


Uc 


7 


4 


0.2424 


1.721 


5 


0.2439 


1.745 


6 


0.2441 


1.749 



Table 5: Linear fits of the sequences {r^, r^ax} for the series (jH]) with assumed fmictional form ((55)) . 
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Figure 13: A comparative plot of estimates of the critical exponent 7 for the Ising model on the CDT 
versus that on a plane triangular lattice (data taken from [36]), obtained via the ratio method. The filled 
shapes refer to the CDT model and the open ones to the plane triangular lattice. The squares indicate fits of 
{rn^i„, r„} sequences, and the circles those of sequences {r„^.^_|_i, r„} with the first ratio eliminated. 

order u^^), even though we cannot claim the resuh to be conclusive because of the limited 
number of terms at our disposal (see Fig. [T3ll . 

What is the picture when we use one of the alternative methods to evaluate the series 
expansions? Results from the Dlog Fade and differential approximants methods are some- 
what inconclusive, most likely because of the small number of estimates we can perform with 
our 6th-order series (see Table [U]). To get a feel for what may be expected at this order, 
we report in Table [7] the corresponding Dlog Fade approximants [(A/" + to the series 

for the plane triangular lattice (see again [36]). The estimates for the critical exponent 7 
in the latter case are clearly closer to the known exact value, whereas for the CDT model, 
with only two values for the diagonal [Af/J\f], it is quite impossible to extrapolate the be- 
haviour of 7. (Note that we are reporting the estimates for the critical point Uc only for 
completeness; they are or course not expected to coincide for the two models.) We have also 
computed the inhomogeneous first-order differential approximants [C/{N' -\- j);N] for the 
CDT series, which was defined just before formula (1451) . Values for the critical exponents 
are not completely off, but there are simply too few of them to make any statement about 
their convergence behaviour (see Table [H]). 



The case of multiple roots can also be considered, but in our case never occurs. 
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Uc 


7 


Uc 


7 
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0.1875 


1.125 


0.2540 


2.064 


0.2513 


2.000 


2 


0.2514 


2.003 


0.2483 


1.8810 


0.2517 


2.008 


3 


0.2519 


2.012 











Table 6: Dlog Pade approximants method applied to the CDT series (|46p . 
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Uc 
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Uc 
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0.2500 


1.500 


0.2666 


1.706 


0.2678 


1.730 


2 


0.2679 


1.732 


0.2670 


1.712 


0.2661 


1.688 


3 


0.2662 


1.692 


0.2667 


1.705 


0.2672 


1.722 



Table 7: Dlog Pade approximants method applied to the plane triangular series in 
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1.396 















1 







Table 8: Inhomogeneous Ist-order differential approximants [C/ {Af + method for the series 
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5.3 The dual CDT lattice series 

The high-temperature series expansion for the Ising susceptibihty on the dual CDT model 
is given by 

XcDTd(M) = 1 + 3m + + 12^3 + 23m^ + (42 + + (78 + + (142 + 

13 1 13 11 

+258m^ + (461 + — + (820 + + (1446 + + (2532 + —)u^^ + 0(n") 

16 8 62 16 

(47) 

Using expressions fl37|) and fl38|) we get the sequence of ratio-method estimates for 7 and Uc 
reported in Table M 



'^max 


Uc 


7 


3 


0.5 


1. 


4 


0.6 


1.6 


5 


0.6147 


1.713 


6 


0.5800 


1.390 


7 


0.5842 


1.436 


8 


0.5782 


1.360 


9 


0.6057 


1.758 


10 


0.6064 


1.769 


11 


0.6093 


1.820 


12 


0.6202 


2.033 



Table 9: Two-point linear extrapolations of the critical point and the critical exponent of the dual series 
(|i7)) obtained using ^ and 

The relative ratio plot is shown in Fig. [HI and shows gentle oscillations around the best linear 
fit of the sequence {r^, ri2}, leading to the estimates 7 = 1.568 and Uc = 0.5956. However, 
these numbers do not tell it all. The strong qualitative differences with the planar honeycomb 
lattice, the appropriate non-random version of the dual CDT lattice, are illustrated in Fig. 
[T5I It shows a comparative plot of the ratios extracted for the two models (the expansion 
coefficients for the honeycomb lattice, known to order u^"^, can be found in [36j). Obviously, 
the infiuence of interfering unphysical singularities present on the regular honeycomb lattice 
is much reduced in the fluctuating dual CDT ensemble. 
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Figure 15: A comparative plot of the ratios of the regular honeycomb lattice (data taken from [35]) versus 
the dual CDT model. 
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On the basis of the rather well-behaved results for the Ising susceptibility on dynamically 
triangulated spacetimes we conjectured in pTS] that "coupling to quantum gravity" may be an 
optimal method to learn about the critical behaviour of a matter or spin system, in the sense 
that the randomness of the underlying geometry eliminates spurious unphysical singularities 
(due to the presence of discrete symmetries in the case of regular lattices), but is not strong 
enough to change the universality class of the matter system. This is a dynamical version of 
similar conjectures made in the 1980s, which advocated the use of fixed random lattices to 
improve the convergence behaviour of lattice gauge theory, say (see, for example, [37|), but 
ultimately did not succeed. 

If it is indeed the case that the singularity structure of the thermodynamic quantities in 
the complex plane is simplified, this could explain that the simple ratio method does indeed 
give the best result, at least at the relatively low orders we have been considering. In the 
case of the Ising model coupled to Euclidean dynamical triangulations, an analogous result 
has already been obtained with regard to the locus of the zeros of the partition function 
in the complex-temperature plane, where unphysical singularities of the flat regular lattices 
have been shown to be largely absent [3H]. However, as already mentioned, the underlying 
geometries are too disordered to serve our present purpose, because their critical matter 
behaviour is changed compared to the flat case. 
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Uc 
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0.5952 


1.554 


0.5574 


0.9738 


0.5686 


1.170 
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0.5691 


1.180 


0.5564 


0.9567 


0.5861 


1.446 
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0.5890 


1.492 


0.5935 


1.560 


0.6359 


2.996 


6 


0.5761 


1.402 











Table 10: Dlog Pade approximants method for the dual CDT series (|47p . 



Returning to the evaluation of the dual CDT series, we have used both the Dlog Pade and 
the differential approximants methods, in addition to the ratio method already described. 
Even with the larger number of terms compared to the original CDT case the results are 
reluctant to show any clear sign of convergence. Our results for the Dlog Pade approximant 
are summarized in Table [TUl The corresponding data for the regular honeycomb lattice [36j . 
which we are including for comparison in Table [TT| seem to show more consistency at the same 
order of approximation. The result from using the inhomogeneous first-order differential 
approximants for the dual CDT lattice (see Table [T2l) is similarly inconclusive. It would 
clearly be desirable to understand in more analytic terms the influence of the fluctuating 
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0.6063 


2.153 
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0.5589 


1.416 


0.5680 


1.531 


0.5676 


1.525 


5 


0.5676 


1.525 


0.5679 


1.530 


0.5589 


1.459 


6 


0.5645 


1.490 


0.5730 


1.618 


0.5709 


1.571 



Table 11: Dlog Pade approximants method for the honeycomb lattice series in [55] , 



geometry of the underlying quantum gravitating lattice on the singularity structure of spin 
models like the one we have been considering, and thus determine which of the approximation 
methods is best suited for extracting their critical behaviour. 

6 Comments on the low-temperature expansion 

Having dealt with the high-temperature expansion, it is natural to ask whether similar ex- 
pansion techniques can be applied in a low-temperature expansion of our coupled models 
of matter and gravity, in the same way this is possible for spin systems on regular lattices. 
The first thing to notice is that at the point of zero temperature around which one expands, 
the spins are all frozen to point into the same direction, and effectively play no role. The 
quantum-gravitational ensemble is therefore again characterized by the probability distribu- 
tion ( ITOi) for the coordination number of pure gravity, as was also the case in the limit of 
infinite temperature. 

In the low-temperature expansion on the regular lattice one starts from a ground state 
with all spins aligned, and then includes perturbations with s overturned spins, obtaining 

m 

In Zn{K, H) = IK + vH + ^[s; r; Gl^V, (48) 

s,r 

in terms of the expansion parameters z = e~'^^ and /i = e'"^^ , where K, H, I and v were 
all defined at the beginning of Sec. 3 above. The part linear in of the (strong) lattice 
constant (see [32] for a definition) for a graph with s vertices and r links embedded in G is 
denoted by [s;r;G], and p is the number of lattice links which are incident on any of the 
vertices of the graph, but do not themselves belong to the embedded graph. On a lattice of 
fixed coordination number q it is easy to prove that p = qs — 2r. Graphs with given values 
of s and r therefore contribute only at a specific order. 

The low-temperature expansion on CDT lattices via graph counting represents a slight 
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Table 12: Inhomogeneous 1st order differential approximants [£/J\f + method for dual CDT series 

(gZl). 
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complication. It is easy to see that already by overturning just a single spin the number of 
links whose interaction changes sign as a result depends on the spin 'sposit ion, so that each 
site will contribute to a different order in the perturbative expansionj^°l In the case at hand, 
this difficulty can be overcome by counting the number of elementary polygons (i.e. faces) 
with given boundary length p on the dual lattice. In the case of several overturned spins we 
can proceed similarly by associating certain patterns on the dual lattice with a given order 
of z. More precisely, to obtain the coefficient at order z^/x* we need to count polygons on the 
dual lattice which have total boundary length p and are constructed out of s faces. Denoting 
the number of such patterns by V{p, s), we can write the free energy per unit volume as 

1 In Z^{K, H) = ^K+^-H + Y, V{p, s)z^^i^ (49) 

s,p 

and the susceptibility at vanishing magnetic field as 

x{z)=Y,^s^^{P^s)zP. (50) 

s,p 

Although the procedure is slightly more involved, and cannot be seen as a computation of 
(averaged) strong lattice constants, the expansion is still doable using the counting techniques 
introduced earlier in this paper. Table [13] gives the results for V{p,s) up to order 10. The 





s = 1 


s = 2 


s = 3 


S = 4 


p 


= 4 


1 

8 








p 


= 5 


1 

8 








p 


= 6 


3 
32 
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3 
128 
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64 


7 
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7 
512 


19 
512 


29 
512 


9 
512 



Table 13: Free energy coefficients V{p, s) as defined in (|49|) for the CDT lattice. 

susceptibility low-temperature series (ISOll turns out to be very irregular and none of the 
methods of analysis considered seems to give a reasonable indication of the critical exponent, 

"'^^By contrast, for the low-temperature expansion on dual CDT lattices, where the coordination number is 
fixed to 3, the expansion parallels that on flat lattices, with the (strong) lattice constants replaced by their 
ensemble averages. 
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at least not from the relatively few terms we have computed. Although this is somewhat 
disappointing, in view of the fact that low-temperature expansions for regular lattices are 
notoriously ill-behaved, it does not come as a total surprise. 

What helps in the evaluation of the low-temperature series for the Ising model on flat 
triangular, square or honeycomb lattices is the fact that their critical temperatures are 
known exactly from duality arguments and thus can be used as an input in so-called biased 
approximants to considerably improve the estimates of critical exponents. Namely, for H = 
one can apply the standard duality transformation \l8l which maps the low-temperature 
expansion ( l48l) of the triangular lattice model to the high-temperature expansion of its 
dual and vice versa. This can be combined with the star-triangle transformation to obtain the 
critical point. The analogous deriviation for the simple square lattice is even easier, because 
it is self-dual. Unfortunately, the star-triangle transformation is not applicable in the CDT 
case, which prevents us from making a similar argument in the coupled Ising-gravity model. 
In order to pursue the analysis of the low-temperature series further, we would therefore 
have to rely on the evaluation of higher-order terms in the expansion. 

7 Conclusions and outlook 

In this article we have given a concrete example of a method which can be used to ex- 
tract physical properties of a strongly coupled model of gravity and matter. We showed 
how the method - estimating critical matter exponents from a series expansion of suitable 
thermodynamic functions of the system - can be adapted successfully from the case of a 
fixed, fiat lattice to that of a fluctuating ensemble of geometries, as is relevant in studies of 
non-perturbative quantum gravity. For the ensemble of two-dimensional causal dynamical 
triangulations, we have formulated an explicit algorithm for counting embedded graphs which 
allows us to do the counting recursively for increasing order. The method can in principle 
be applied to other matter and spin systems which admit a similar diagrammatic expansion 
in terms of weak or strong lattice constants around infinite or zero temperature, and where 
one has sufficient information about the probability distribution of local geometric lattice 
properties like the vertex coordination number. Even in cases where these lattice properties 
are not available exphcitly, the lattice constants could still be extracted from simulating the 
pure gravity ensemble. 

As a potential spin-off, we noticed that the explicit expansion for the susceptibility in 
vanishing external field of the Ising model on CDT (up to order 6) or dual CDT lattices (up 
to order 12) indicates a more regular behaviour of this function in the complex-temperature 
plane than that of the Ising model on the corresponding fixed triangular or honeycomb lattice. 
This may also explain why the simple ratio method, applied to CDT lattice results from only 
six orders, gives an excellent approximation to the Onsager susceptibility exponent 7 = 1.75. 
As we saw in some detail, other approximation methods, namely, Dlog Fade and differential 
approximants, do not produce results of a similar quality. We believe that for the case of the 
dual CDT lattice, the series computed is simply too short to yield reliable estimates with any 
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of the approximation methods (even at order 12, the number of terms contributing is much 
smaller than the number of terms contributing at order 6 on the original CDT lattices). At 
any rate, if one aims to make an argument of improved convergence of matter behaviour on 
fluctuating lattices (see also [12]), it is plausible that this will be achieved optimally with 
the triangulated CDT geometries whose vertex coordination number can vary dynamically 
from 4 all the way to infinity, rather than with the dual tesselations which have a fixed, low 
coordination number of 3. To investigate the convergence issue in more detail will require 
going to higher order than 6 on CDT lattices, which is not really feasible 'by hand', as we 
have been doing so far, but will require the setting up of a computer algorithm to perform 
(at least part of) the graph counting. 

Our work should be seen as contributing to the study of non-perturbative systems of 
quantum gravity coupled to matter, which is only just beginning. One challenging task will 
be to establish computable criteria characterizing and quantifying the influence of geometry 
on matter and vice versa, in the physically relevant case of four spacetime dimensions, 
something about which we currently know close to nothing. Obvious quantities of interest 
are critical exponents pertaining to geometry (like, for instance, the Hausdorff and spectral 
dimensions of spacetime already measured for the ground state of four-dimensional CDT 
[7j), and critical matter exponents, like that of the susceptibility investigated in the present 
work. One would like to have a classification of possible universality classes of gravity- 
matter models as a function of the characteristics of the ensemble of quantum-fluctuating 
geometries. A necessary condition for viable quantum gravity models is that at sufficiently 
large distances and for sufficiently weak matter fields the correct classical limits for these 
critical parameters must be recovered. 

Two-dimensional toy models like the one we have been considering can serve as a blueprint 
for what phenomena one might expect to find. In two dimensions, there have been several 
investigations of the effect of random geometric disorder on the critical properties of matter 
systems, and attempts to formulate general criteria for when a particular type of disorder is 
relevant, i.e. will lead to a matter behaviour different from that on fixed, regular lattices. 
Good examples are the Harris [10] and the Harris-Luck [T3] criteria, which tie the relevance 
of random disorder to the value of the specific-heat exponent a and to correlations among 
the disorder degrees of freedom. However, not all models fit the predictions, and open 
problems remain (see [T6] and references therein). Causal dynamical triangulations coupled 
to matter add another class of models inspired by quantum gravity, whose randomness lies in 
between that of Poissonian Voronoi-Delaunay triangulations and the highly fractal Euclidean 
dynamical triangulations, and over which one has some analytic control. It will be interesting 
to see to what degree the robustness of the matter behaviour with respect to the geometric 
fiuctuations observed so far [TOl [TTl 112] will persist for different types of spin and matter 
systems and in higher dimensions. 
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